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This is a historical note. Bethe Ansatz solvable models are considered, like XXZ Heisenberg 
anti-ferromagnet and Bose gas with delta interaction. Periodic boundary conditions lead to Bethe 
equation. The square of the norm of Bethe wave function is equal to a determinant of linearized 
system of Bethe equations (determinant of matrix of second derivatives of Yang action). The proof 
was first published in Communications in Mathematical Physics, vol 86, page 391 in 1982. Also 
domain wall boundary conditions for 6 vertex model were discovered in the same paper [see 
Appendix D]. These play an important role for algebraic combinatorics: alternating sign matrices, 
domino tiling and plane partition. 



INTRODUCTION 

Many two dimensional models have been solved by means of the Bethe Ansatz, see for example [H, 0, S 0, S| (30| - 
Quantum inverse scattering method (QISM) [fj, |7|, l26[ discloses algebraic nature of these solutions. Michael Gaudin 
studied norms of Bethe wave functions for the quantum nonlinear Schroedinger equation 31 1 Q and suggested a 
remarkable conjecture that the norm of the eigenfunction is equal to a Jacobian. In [25| these formulae are proved 
and a more general result is obtained. The norms are calculated for any exactly solvable models with the R matrix 
either of the XXX model or of the XXZ Heisenberg models. This note gives an idea about an approach of the 
paper. 



NOTATIONS 



First of all let us remind the reader of some notations of the QISM. Eigenfunctions of the Hamiltonian of the 
physical system are constructed by means of the monodromy matrix of an auxiliary linear problem T(X). In our case 
T(A) is a 2 x 2 matrix, its matrix elements being quantum operators, which depend on the spectral parameter A: 

T(A) - C(A) D(X) ) 
Commutation relations of these entries are given by the formula: 



R(X, n) (T(A) <g> I) (I ® T(/i)) = (/ <8> T(n)) (T(A) <g> /) R(X, (i) (1) 

Here / is the unit 2x2matrix, i?(A, /i) is a 4x4 matrix with c-number elements. Another way to write this equation is: 

R(X, n) a ,pT{\) a T{ji) p = T{fi) p T{X) a R{X, n) a , (2) 

This means hat T(X) a is a 2 x 2 matrix which acts on the 2 dimensional space with index a and T(/j,)p is a matrix 
which acts on space with index (3. We assume that a ^ (3. The matrix R(X, )J.) a ,(3, acts in the tensor product of these 
two spaces. We shall deal with R matrices of the following form: 



R(X,n) = 



//(A*, A) 
1 

g(fx,X) 
\ 



\ 





g(fx,X) 

1 

o /(m,a)7 



For models of the XXZ type: 



/(a,aO = 



For models of XXX type: 



/(A, A*) 



sinh(A — fj, + 2ir/) 
sinh(A — fi) 

(A — fi + 2ik) 



ff(A,/i) = 



i sin(2?7) 
sinh(A — /i) 



(A - M) 



(a-m) 



(3) 



(4) 
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Here A and \i are spectral parameters and r\ and k are coupling constants. Note that XXX is a limiting case of XXZ. 
We can also write the R matrix using notations of formula ([2]). 



R(X,») a , = cosfo) ^fe^ffl I a Ip + i Ml) C °^^ oi4 (5) 

(6) 

Here a are the standard Pauli matrices and 2cr ± = a 1 ± a 2 . Let us write down some commutation relations of (TT|) 
explicitly 

[B(X),B(fj,)} = 0, [C(A),C(/i)]=0, [A(A) + £>(A),A( M )+ J D( M )] = (7) 
A( M )£(A) - /(/*, A)S(A)A( M ) + g(X, ii)B(p)A(X) (8) 
D(Ai)B(A) - /(A, Li)B(\)D(p) + g((i, X)B^)D(X) (9) 
[C(A C ), i?(A b )] = 5 (A C , X b ){A(X c )D{X b ) - A(X b )D(X c )} (10) 

Other commutation relations can be found in the book [26[. Pseudo- vacuum |0 > and dual pseudo- vacuum < 0| are 
important: 

C(A)|0>=0, A(A)|0 >= o(A)|0 >, D(A)|0 >= d(A)|0 > (11) 
< 015(A) = 0, < 0|A(A) =< 0|a(A), < 0|D(A) =< 0|d(A) (12) 

Here a(A) and d(X) are complex valued functions. The space in which the operators A(X), B(X), C((A) and D(X) act 
is constructed in [26|]. The norms and scalar products in question are functionals of these a(X) and d(X). We shall 
vary a(A) and d(X) and study the dependence of scalar products on these functional arguments. The Hamiltonian of 
the physical system in question is expressed in terms of the transfer matrix t(/i) = A(X) + D(X). The Hamiltonian 
and the transfer matrix have common eigenfunctions which are constructed as follows. Put 

i> N ({X j }) = B(X 1 )...B(X N )\0> (13) 

and suppose that Xj satisfy the system of the transcendental equations (TE) : 

a(X n ) tt /(A n , Xj) . , . 

^n/(A-xo =i here 3+n (14) 

Then ipN({Xj}) is an eigenfunction of t(A) with the eigenvalue: 

JV N 

B(ji) = a(n) J] ffa, Aj ) + d{n) J] f(Xj , M ) (15) 

j=l 3=1 

Here N is called the number of particles. Note that 



ip N ({X j })=<0\C(X 1 )...C(X N ) (16) 

is a dual eig enfunction for t(fi) with the same eigenvalue. Pauli principle was proved in the original publication [25| . 
see also 26|. We take all Xj to be different. Finally let us present two remarks. First of all new variables 

* k = ^^) + \^nx^xTy k ^ 1 -- N (17) 

are convenient. For example (|14[) can be written as 

(/) k = mod27r (18) 
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EXPRESSION FOR THE NORM IN THE XXX TYPE MODELS 

For models with an R matrix ([4]) the scalar product of an eigenfunction and dual eigcnfunction is equal to 

< 0|C(Ai) • . . C(X N )B(X 1 ) . . . B(X N )\0 >= k n (]jf =1 a(Xj)d(Xjj) < > (19) 



(nf =i nf =lj¥fc /(A„A fe )) det N (|£ 



Here the cj> k are the variables (fTT|) and the set of the Xj is a solution of the system 1(15]). (fT4")) .TIie derivatives can be 
written in the explicit form: 

f¥fe = c (, d a(A fc ) A 2k \ 2k 

d\j k ' 3 \ dx k d(x k ) + ^ (x k - x p y + K*J (A fe - Xj) 2 + « 2 1 j 



The formula (|19[) is useful for calculation of norms if the set of Xj invariant under complex conjugation {Xj} — {Xj} 
and -B' (A) = ±C(A). Examples of applications can be useful. Let consider nonlinear Schroedinger equation, it is also 
known as Bose gas with delta interaction. It has a Hamiltonian 

H = J dx (d x ip^d x ip + nift ipip) > 

[l>(x),tf(y) = 5(x-y)] (21) 

The monodromy matrix has the following property: 

a x T^(X)a l = T(A), St (A) = (7(A) (22) 

This is XXX case. The vacuum eigenvalues are equal to: 

o(A) = exp(-iLA/2), d(X) = exp(iLA/2) 

This proves Gaudin conjecture Q. This determinant formula for norm also applicable for Heisenberg XXX spin chain. 

EXPRESSION FOR THE NORM IN THE XXZ TYPE MODELS 

For solvable models with XXZ R matrix ^ the scalar product of the eigenfunction and dual eigcnfunction is equal 
< 0|C(Ai) . . . C(X N )B(Xi) . . . B(X N )\0 >= (sm2 v ) N (j[f =1 a(X 3 )d(X^ < 0|0 > (23) 
(uUUtx,j^f(^h)) delar(Sfc) 
Here Xj has to satisfy equation (1 14|) . The Jacobi matrix can be written down in the explicit form: 

*b k r ( .d\na(X k )/d{X k ) A sin(4?7) \ - sin(4r?) 



<JVk _ e . u iuayA k )/u{A k ) \ - am^tr/j 

a\— k,j \l „ x + 2^ cinlnQ, _ \ l OinUinlna,. - \ . - 9™^ sir 



dXj ' \ dX k _ sinh(Afc — A p + 2irj) sinh(Afc — X p — 2irj) I sinh(Afc — Xj + 2ir\) sinh(Afc — Xj — 2ir\) 

1 (24) 
The whole paper [25| is devoted to the derivation of this formula. Let us calculate the norms for X XZ model. The 

Hamiltonian of the model is 



M 



H = J2 a ^l+i + 44+1 + cos(2r])(a 3 k al +1 - 1) (25) 
fe=i 

The model was imbedded into QISM in [7]. This monodromy matrix has the following property at real rj: 

a 2 T+(A> 2 = T(A), 5+(A) = -(7(A) (26) 
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Pseudo-vacuum is the ferromagnetic state |0 >= J^J -_j Tj = < 0|. The vacuum eigenvalues are: 

<x(A) = sinh M (A — irj), d(X) = sinh (A + irj) 
In order to write down the square of the norm it is convenient to introduce: 

sin(2r;) 



sinh(A — irj) sinh(A + in) 



The formula for the square of the norm is: 



< 0|flt(A JV ) . . . St(A x )S(Ai) . . . B(X N )\0 >= sin N (2r)) ( sinh AJ (A-,- - irj) sinh 717 (Aj + irj) 



(rifc>j=i 



sinh(Aj —\k—2iri) sinh(Aj — Afc+227/) 
sinh 2 (Aj— A fc ) 



(27) 



) dctjv [S kj (M X {X k , v) ~ Eili x(A fc - A z , 2??)) + x(A fc - Xj,2r)) 



The formula was presented in [9( and verifies for N = 2 and N = 3. The proof for arbitrary N was first published in 
[2BI I . The formula (|2"3"|) also describes norms in Sine-Gordon and lattice Sine-Gordon 18, 



THE IDEA OF THE PROOF 

In order to prove determinant formula for the norm of Bethe wave function let us introduce an object: 

< 0|C(Ai) . . . C(Aiv)S(Ai) . . . S(Ajy)|0 > 



( S in2^ [nu'&Wi)) <oio> (n- v =inLi,^ fc /(^,A fc ) 

We have to prove that 



Al..A 



IA1 ... A/vl = det 



n \d\ k 



(28) 



(29) 



We will assume that < 1 >= 1. The author of [25( proved that a(A) and d(X) can be considered an arbitrary 
functions. So we can consider the variables 

a{X p ) 

Ap - l dx p [n d(x p ) 



as independent of Aj , see [25( . The following theorem is proved in [25( . 

In order to prove i29\) in is enough to prove the following five properties of \X% . . . X^\: 

• It is invariant under simultaneous replacement 

Xj «-> Xk, Xj «-> Xk 

• It is leaner function of X\ 

\X x ...X N \ = UiX x + Vi 

• The coefficient U\ is 

Ui = \x 2 . . . x N \ modl fi ed 



(30) 
(31) 
(32) 



The right hand side is given by formula (28\) with X± removed and functions a(A) and d(X) replaced by: 

a modified {x) = a ( A)/ ( A) Al)j d modified {x) = d{X )f(X 1 , A) (33) 

• It vanish if all X p = 

|A x ...Aiv| = 0, if all X p = 0, at p = l,...,N (34) 

• For N = 1 

I Ai| =X X 



The paper [25| proves that the right hand side of ([28)1 has all five properties listed above. The proof is reduced to 
analysis of six vertex model with domain wall boundary conditions. 
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SIX VERTEX MODEL 



The six vertex model is an important 'counterexample' of statistical mechanics: the bulk free energy depends on 
the boundary conditions even in thermodynamic limit, see [ill ]. 
Let us start the formal presention: L operator in site number k is 

/ sinh(A - v k ~ irial) -ia^ sin(2?7) 
k{ k> \ -ia+ sin(277) smh(A - v k + iwl) 

It also can written as: 

Lk(X — Vk) = cos(?7) sinh(A — v k ) ~ i singer 3 cosh(A — v^) — i sm(2i])(a + a^ + a~<J k ) (35) 
The monodromy matrix: 

T(A) = L M (A - v M ) ...Lt(X- Vl ) (36) 
obeys the commutation relations |T]) with R matrix ([3]). Pseudo- vacuum is the ferromagnetic state 

M 

io >= n ti=< oi 

The vacuum eigenvalues are now equal to 

M M 

a(A) = J^J sinh(A — fj — iff), d{X) = J^J sinh(A — Vj + ii]) 
3=1 3=1 

Such an inhomogeneous generalization was used for example in [20l 21 , 22[ ■ Let us consider a special state 

B(X 1 )...B(X M )\0> (37) 

Here the number of the B(X) is equal to the number of the sites in the lattice N — M. All spins are looking down in 
this state: 



Here Zn is a complex number and 



B(X 1 )...B(X M )\0>= Z N Q (38) 



= 1] ii ( 39 ) 

3=1 



The definition of Zn is 



Z N = CIB(X 1 )...B(X M )\0> 

The paper [25j proves that Zn is the partition function of six vertex model with domain wall boundary conditions, 
see Appendix D. Actually these boundary conditions were introduced in this paper, including the name. Explicit 
description of six vertex model with domain wall boundary conditions also can be found in [2II ] . 

The following recursion relations were discovered in [25[: If v\ = Ai + if] then Zjq reduces to Zm-i with Ai and v\ 
removed. 

Z N ({X a }, {fjJOU-A^iT) = -«sin(2r?) Y\_ sinh(Ai - v k - ir]) J ]^[ sinh(A Q - z/ x - £77) Zat-i({A q ^i}, {vj^x}) 

\k=2 ) \a=2 ) 

The derivation of this recursion relations also can be found in section 2 of the paper [2S| . 



6 



CONCLUSION 

There was a lot of progress since the paper 12511 . The partition function of six vertex model with domain wall 
boundary conditions has a lot of applications [32i | and generalizations, see for example [29j]. Still there are open 
problems. For example the determinant formula for Bethe wave function in Hubbard still not proven, see [2^ ]. 
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